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ABSTRACT 

Whether or not supernovae contribute significantly to the overall dust budget is a contro- 
versial subject. Submillimetre (submm) observations, sensitive to cold dust, have shown an 
excess at 450 and 850 /im in young remnants Cassiopeia A (Cas A) and Kepler. Some of 
the submm emission from Cas A has been shown to be contaminated by unrelated material 
along the line of sight. In this paper we explore the emission from material towards Kepler 
using submm continuum imaging and spectroscopic observations of atomic and molecular 
gas, via Hi, 12 CO (J=2-l) and 13 CO (J=2-l). We detect weak CO emission (peak T£ = 
0.2-1 K, l-2kms _1 FWHM) from diffuse, optically thin gas at the locations of some of the 
submm clumps. The contribution to the submm emission from foreground molecular and 
atomic clouds is negligible. The revised dust mass for Kepler's remnant is O.1-1.2M , about 
half of the quoted values in the original study by Morgan et al. (2003), but still sufficient to 
explain the origin of dust at high redshifts. 
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1 INTRODUCTION 

The conditions following a supernova (SN) explosion are thought 
to be conducive to the formation of dust: the abundances of heavy 
elements are high, as is the density; temperatures drop rapidly in 
the expanding ejecta, quickly reaching levels allowing the subli- 
mation of grain materials. Theoretical estimates predict that type-ll 
SNe should produce a significant quantity of dust, approximately 
0.3 — 5 M© per star, depending on the metallicity, stellar mass and 
energy of the explosion (e.g. Todini & Ferrara 2001; Nozawa et al. 
2003; Schneider, Ferrara & Salvaterra 2004). Some circumstantial 
evidence also leads us to believe that SNe should be an important 
source of dust: first, without SNe there is a dust budget crisis in 
the Galactic ISM. The dust produced in cool stellar atmospheres 
of intermediate-mass stars, combined with current predictions for 
how much dust is destroyed in shocks, yields far less dust than is 
observed (Jones et al. 1994) in the interstellar medium (ISM). Ei- 
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ther another source of dust exists, or dust destruction cannot be 
as efficient as is widely believed. Second, without SNe and their 
massive precursors as significant sources of dust, it is difficult to 
explain the immense dust masses found in submillimetre- selected 
galaxies and quasars at high redshift (e.g. Smail et al. 1997; Isaak 
et al. 2002; Eales et al. 2003). There is not sufficient time for dust 
to form in such large quantities from evolved stars alone (Morgan 
& Edmunds 2003; Dwek et al. 2007 and references therein). 

The signature of warm, freshly-formed dust in Cas A was 
seen in spectroscopic data by Rho et al. (2008), with reported 
dust masses in the range O.O2-O.O5-M . Smaller fractions of warm 
dust have been reported in SN2003gd and Kepler (Sugerman et al. 
2006; cf. Meikle et al. 2007; Blair et al. 2007). Spitzer observations 
of SNRs in the Magellanic Clouds are also consistent with small 
amounts of dust (e.g. Borkowski et al. 2006; Williams et al. 2006). 
However, Spitzer is not sensitive to the presence of very cold dust, 
which peaks at wavelengths longer than 160 /mi. To address the 
question of whether large quantities of dust are present, we require 
observations at longer wavelengths in the submm. 

The Submillimetre Common User Bolometer Array (SCUBA 
- Holland et al. 1999) on the James Clerk Maxwell Telescope 
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(JCMT0) was used to observe the young Galactic SN remnant 
(SNR) Cas A (Dunne et al. 2003 - hereafter D03), with large ex- 
cesses of submm emission detected over and above the extrapolated 
synchrotron components. This was confirmed by ARCHEOPS 
(Desert et al. 2008). Because of the high spatial correlation with 
the X-ray and radio emission, this was interpreted as emission 
from cold dust associated with the remnant. However, some of the 
submm emission comes from molecular clouds along the line of 
sight (Krause et al. 2004; Wilson & Batrla 2005). The high degree 
of submm polarisation suggests that a significant fraction of dust 
does originate within the remnant (~ 1 Mq, Dunne et al. 2009). We 
now turn to the only other remnant with a reported excess of submm 
emission over the extrapolated synchrotron: Kepler. We originally 
interpreted our SCUBA data as evidence for 0.3-3 M of dust as- 
sociated with the remnant (Morgan et al. 2003 - hereafter M03; 
Gomez et al. 2007), an order of magnitude higher than those pre- 
dicted from 160-/xm Spitzer data of Kepler (Blair et al. 2007). 

Kepler's SN has a shell-like structure, ~3 arcmin in diameter. 
Estimates of its distance, using H I absorption features range from 
3.9 to 6kpc (Reynoso & Goss 1999 - hereafter RG99; Sankrit et 
al. 2008). Its classification has been controversial (Blair et al. 2007; 
Reynolds et al. 2007; Sankrit et al. 2008) with evidence pointing to- 
wards either a type-la - the thermonuclear explosion of a low-mass 
accreting star in a binary system - or a type-lb - the core collapse 
of a massive star. Reynolds et al. proposed that Kepler's SN was 
the result of a thermonuclear explosion in a single 8-M© star, af- 
ter roughly 50 Myr of evolution, which would make it pertinent to 
the issue of the dust budget in the early Universe. Of course, un- 
derstanding dust formation in SNRs is important regardless of the 
explosion mechanism: dust formation following type-la SNe would 
indicate that type-ll SNe would also be likely dust producers (Clay- 
ton et al. 1997; Travaglio et al. 1999). Here, we revisit the submm 
data and ask how much of the submm emission can be associated 
with the remnant and how much with material along the line of 
sight. 

In this paper, we recalculate the contribution of synchrotron 
radiation to the submm emission in Kepler using a radio spectral 
index map (kindly provided by T. DeLaney). We also present the 
first high-resolution molecular-line map towards Kepler's SNR and 
compare this with archival Hi data (RG99). In fl2.1l we present de- 
tails of the observations and data reduction; in ^2.2l we investigate 
the possibility of contamination by line-of-sight gas clouds. In ^3] 
we compare the distribution of the dust and gas. In ^3.11 we esti- 
mate the dust mass in Kepler with results summarised in Sim- 
ulations of the effects of the SCUBA chop throw are discussed in 
Appendi?j4] 



2 OBSERVATIONS AND ANALYSIS 

2.1 Submillimetre continuum observations 

The reduction of SCUBA data of Kepler's remnant at 450 and 
850 /xm, using standard routines in the SURF software package 
(Sandell et al. 2001), was described briefly by M03. The obser- 
vations were carried out over five different nights during 2001- 

03 using the jiggle-map mode. The array was chopped to remove 
sky emission; we chose reference positions using a map of the 

1 The JCMT is operated by the Joint Astronomy Centre on behalf of the 
UK's Science and Technology Facilities Council, the Netherlands Organisa- 
tion for Scientific Research, and the National Research Council of Canada. 



radio emission as a guide, with a chop throw of 180arcsec (see 
Fig. lAlb . The potential hazards of this observing mode are dis- 
cussed and simulated in In M03 we presented the synchrotron- 
subtracted 850-/xm image, obtained by scaling the 5-GHz Very 
Large Array (VL43) map of Kepler using a constant spectral in- 
dex, S u oc , that being the mean value reported by DeLaney 
et al. (2002). We take a different approach here, instead using the 
spectral index map of DeLaney et al. - their Figure 4 - to produce 
a more accurate estimate of the synchrotron contribution pixel-by- 
pixel. The flatter spectral index in the north of the remnant means 
that more flux has been subtracted from the submm images than 
in the analysis of M03. The final synchrotron- subtracted submm 
flux densities are 0.4 ± 0.1 and 2.8 ± 0.7 Jy at 850 and 450 /xm, 
respectively. 

The signal-to-noise maps of cold dust at 450 and 850 /xm are 
shown in Fig. [TJ with the synchrotron contribution subtracted us- 
ing the spectral index image. The noise images were created by 
randomly generating 1,000 artificial images (Eales et al. 2000). 
The noise is substantially higher near the edges of the map be- 
cause these regions received significantly less integration time. At 
450 /xm, the centre of the map is also noisy (the array footprint is 
smaller at 450 /xm than at 850 /xm). The peak signal-to-noise val- 
ues in these maps are 5a at 450 /xm and la at 850 /xm. There are 
regions of high flux near the edges of the map (particularly in the 
south-east) but these are also regions of high noise so our confi- 
dence in these features is low. Conversly, there are regions of low 
flux which are not particularly noisy (e.g. in the southern 'ring' re- 
gion of the remnant); our confidence in these features is commen- 
surately higher. 

The location of the submm peaks A-J, are marked on the 1.4- 
GHz map in Fig. [2 This image was constructed using data from 
the VLA archive (see ^2.2.31 ). Dust clumps A, B, F, I and G all 
fall within the radius of the shockfront. Most of the emission from 
cloud E is within the shock also whereas C, D, H, and J lie beyond 
the X-ray and radio boundaries. Cloud E is located at the position 
of one of the 'ears' in the radio, a region associated with the ejecta 
but beyond the almost circular shock front. Around 30-40 per cent 
of the submm flux lies outside the shock front (as defined by the 
radio observations at lOOarcsec). 

Possible explanations for the submm emission are: (i) dust 
produced by the supernova remnant or progenitor star; (ii) inter- 
stellar material along the line of sight and/or (iii) spurious structure 
in the SCUBA map - perhaps an artefact of the observing or data 
processing techniques. These are possibilities we will explore in 
the remainder of this paper. 

2.2 Exploring the possibility of line-of-sight contamination 

Does Kepler's SNR have significant interstellar material along the 
line of sight that may be contributing to the measured submm 
fluxes? Kepler's SNR is approximately 600 pc out of the Galactic 
Plane (/ ~ 4.5°, b ~ 6.8°) with extinction and 100-/xm back- 
ground lower than measured at Cas A by a factor of three (IRAS 
IRSKY maps - Arendt 1989). A latitude-velocity map of integrated 
12 CO ( J=l-0) (kindly provided by T. Dame) at the location of Ke- 
pler shows that the velocity range of clouds integrated over longi- 
tude is confined to — 5 < v < 5kms _1 (Fig. [3^). Fig. [3]} shows 

2 The VLA is operated by the National Radio Astronomy Observatory, 
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operative agreement by Associated Universities, Inc. 



Foreground contribution of dust towards Kepler 3 






17:31:00 55 50 45 40 35 30:30 25 



RA (J2000) 

Figure 2. 1.4-GHz VLA image of Kepler's SNR. The rms noise is 
0.6 mJy beam -1 . Submm clumps are labelled A-J with radius, 23arcsec. 
The green contour indicates the boundary of the SCUBA observations and 
the black box indicates the extent of the molecular observations in this work 

(§ED- 

the integrated CO emission near Kepler's SNR (marked by the ra- 
dio contour in the centre) from the processed Galactic CO survey 
(Dame et al. 2001). The remnant is within 3° of the p Oph cloud 
complex (distance, 165 pc) which has clouds with line widths (Av) 
up to 100 kms -1 . 

2.2. 1 12 CO( J =2-1 ) observations: searching for molecular 
structures 

In order to quantify the possible contribution from foreground 
molecular material, we observed Kepler's SNR in the 12 CO( J=2- 
1) line using the A3 receiver on the JCMT with the Auto- 
Correlation Spectrometer Imaging System (ACSIS) in double- 
sideband mode (DSB) (Table [TJ. The bandwith was 1.8 GHz with 



1,904 channels and a velocity coverage of 2,000 km s _1 . We raster- 
mapped over a region 8x8 arcmin 2 with grid points separated by 
7 arcsec to avoid smearing in the scan direction (the beam at this 
wavelength is 20.8 arcsec FWHM). The integration time was 5 s at 
each point. In this intial dataset, the lines were unresolved. Fur- 
ther high-resolution service observations were obtained to resolve 
the lines at the location of two CO peaks. The high-resolution data 
were taken with the same receiver set-up, but with a bandwidth of 
250 MHz, 8,192 channels and a velocity coverage of 300kms _1 . 
Spectral resolutions and positions are given in Table Q] The spectra 
were baseline- subtracted with a third-order polynomial and scaled 
to main beam (Tmb) temperatures by dividing the antenna temper- 
atures (TX) by the telescope efficiency ?7mb ~ 0.69. 

We detect faint, narrow emission corresponding to the radial 
velocities of the CO and H I emission in the nearby p Ophiuchus 
cloud complex between — 8 < v < 20kms _1 (de Geus 1992). 
Channel maps showing the total CO brightness temperature over 
the range —7.6 < v < 20.3kms _1 with velocity intervals of 
1.9 km s -1 are shown in Fig. [4] This range encompasses all of the 
CO lines detected over the entire velocity coverage. Lines are ob- 
served at velocities of —1.4, 3.8 and 11.5kms _1 with linewidths 
(FWHM) Av < 0.8 km s -1 , indicating cold, faint clouds. These 
lines are narrower than expected for giant molecular clouds, typ- 
ical interstellar material (GMCs; Solomon, Sanders & Scoville 
1979, Av ^4-7 km s -1 ), so-called dark clouds and high-Galactic- 
latitude clouds (Lada et al. 2003). The high-resolution spectra at 
locations listed in Table Q] were fitted with Gaussian profiles using 
the SPLAT-VO package (Draper & Taylor 2006) with CO intensi- 
ties ranging from 0.5-1.4 K km s -1 . It is difficult to ascertain dis- 
tances to these structures since the Galactic rotation models break 
down towards the Galactic Centre. Without an estimate of distance 
to the CO structures, the observed size- velocity dispersion relation- 
ship for molecular clouds in the Milky Way (Tsuboi & Miyazaki 
1998; Oka et al. 1998) suggests that the CO clouds are smaller than 
1 pc across. 

We see no evidence for shocked molecular gas which would 
broaden the lines to average linewidths Av>15kms _1 , as seen in 
the SNRs W28, W44 and W51C (De Noyer 1979; Junkes et al. 
1992; Koo & Moon 1997; Seta et al. 1998; Reach, Rho & Jar- 
rett 2005). If an interaction is present, we would expect to ob- 
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Figure 3. (a) Latitude-velocity map at a longitude of 4.5° from the CO Milky Way Survey (e.g. Dame et al. 2001). Contours show T£ ranging from 0.3 to 
4.5 K in steps of 0.6 K. A line at a constant latitude of 6.8° is included to show the location of Kepler's SNR. The greyscale goes from to 8.5 K, as shown at 
the right of Fig. 3b. (b) Integrated CO image of the molecular clouds (part of the Ophiuchus complex at 165 pc) in the vicinity of Kepler's SNR, taken from the 
Galactic CO data available at http : //cfa-www . harvard . edu/mmw/MilkyWaynMolClouds . html The resolution is 8 arcmin and the integration 
spans over —300 < v < 300 kms -1 . The contours are at constant T£, ranging from 1 to 9 K in steps of 1 K. The greyscale is indicated to the right of the 
image. An outer radio continuum contour of Kepler's SNR with location indicated by the arrow. 
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Figure 4. Channel maps of Tj£ across the 12 CO (J=2-l) transition ranging from —10 to 24.6 kms -1 with velocity width, Av = 1.9kms -1 . The images 
have been smoothed with a 21-arcsec gaussian. The central velocity is labeled on each grid and the negative greyscale ranges from —0.41 to 1.2 K kms -1 . 
Three clouds are detected: to the south of the remnant at velocities 3-5 kms -1 , to the north at velocities 9-11 kms -1 and a small cloud in the north at 
velocities 0-2 kms -1 . 
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Figure 5. Grid of 12 CO (J=2-l) spectra towards Kepler's SNR from —5 < v < +15kms _1 across the entire area. The black contour defines the shockfront 
as estimated from radio observations (Fig. [3. The faint, narrow-line emission from the two molecular clouds in the south and in the north are clearly visible 
outside the area encompassed by the supernova remnant. 



serve CO emission peaking at the location of the shock front (e.g. 
Wilner et al. 1998) and a velocity jump across the outer shock 
front. Fig. \5\ shows the grid of spectra across the entire map, 
smoothed to a velocity resolution of 3kms _1 and restricted to 
— 5 < v < +15kms _1 for clarity. The location of the super- 
nova shockfront (defined by the outer radio contour) is overlaid. 
The molecular clouds detected here lie outside the shock front and 
line-broadening is not observed in the profiles. We also investigated 



the velocity gradients across the clouds by integrating in R.A. and 
declination (Fig. [6]). The CO clouds at 4 and llkms -1 line up 
at a common R.A., spanning 280arcsec with a velocity range of 
~2kms _1 . The clouds span a similar range in Dec. and velocity 
in the RA. velocity-position map; the lack of velocity gradient in 
these images confirms there is no interaction. 
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Line 


R.A. 


Dec. 


t (mins.) 


Date (2007) 


Beam (hpbw) 


Spectral resolution 


Pos. switch to: 


12 CO (J=2-l) 


17:30:42.0 


-21:29:35.4 


1320 


March 08 - April 08 


20.8" 


977 kHz, 1.27 km s -1 


-426.0", +72.0" 


12 CO (J=2-l) 


17:30:45.8 


-21:33:10.9 


57 


June 04 


20.8" 


31kHz, 0.08 kms- 1 


+163.8", +374.1" 


12 CO (J=2-l) 


17:30:46.0 


-21:27:23.7 


167 


June 04 


20.8" 


31kHz, 0.08 kms -1 


+161", +721.3" 


13 CO (J=2-l) 


17:30:45.8 


-21:33:10.9 


167 


June 04, 08 


21.3" 


31kHz, 0.08 kms" 1 


+163.8", +374.1" 



Table 1. Summary of the CO observations used in this work. 
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Figure 6. High-resolution 12 CO (J=2-l) spatial-velocity maps of the re- 
gion surrounding Kepler's SNR. Top: R.A. -velocity map with contours 
starting at 0.04 K and incrementing by 0.03 K. Bottom: Dec. -velocity map 
with contours starting 0.07 K and incrementing by 0.03 K. 



2.2.2 16 CO(J=2-l ) observations 

We observed the brightest 12 CO(J=2-l) cloud, at 11 kms -1 , in 
13 CO ( J=2-l) and detected no emission; the high-resolution spec- 
trum is shown in Fig. [7] The 3 a upper limit can be estimated us- 
ing ico < 3ctTmb V Av 5v, where Sv is the velocity width of the 
channel and Av is the velocity width of the expected line. We 
find an upper limit of /( 13 CO) < 0.08 Kkms" 1 , for a width 
of lkrns" 1 (the maximum width detected at the same location 
in 12 CO). The lower limit on the brightness ratio R is therefore 
7( 12 CO)/^( 13 CO) >17. This is larger than those usually found 
in GMCs where R - 5 (Solomon et al. 1979). 

The optical depth, tJ 3 , is estimated using Eq.Q] 
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Figure 7. 12 CO(J=2-l) high-resolution spectra (black) observed at the 
location listed in Table[T]over the velocity range — 14 < v < 23 kms -1 . 
The corresponding 13 CO(J=2-l) spectra are shown in purple. The peak 
CO emission is <1 with velocity widths Av <1 kms -1 . 



perature of the gas T ex , we assume a nominal value of 15 K (similar 
to the dust temperatures implied by the 450-/850-/im flux ratios, 
M03). The upper limit on the optical depth is rj 3 < 8 x 10~ 3 for 
T ex =15 K, assuming that T ex is the same for the 12 CO and 13 CO 
molecules. The non-detection of 13 CO indicates that the 12 CO 
emission is optically thin and confirms that we are observing faint 
emission from a diffuse molecular cloud. Under the assumption that 
the excitation temperature is the same for both molecules and emis- 
sion fills the beam, the 3-cr upper limit on the optical depth of the 
12 CO emission is rj 2 < 0.15 (Eq.EJ). 



2 CO), 



exp 



T mb ( 13 CO). 



exp (- 



(2) 



where To=h v jk. Without further knowledge of the excitation tem- 



The beam-averaged total column density, N, can be estimated 
using the standard assumptions of radiative transfer in local ther- 
modynamic equilibrium (LTE) for an optically thin gas (Thompson 
& MacDonald 1999): 

3k fTuBdv / E u \ 

N ~Q~l\~Q~ll Q(^ex)eX P y— - , (3) 

where v is the line frequency, S is the line strength, fx is the perma- 
nent electric dipole moment (in e.s.u), gi and #k are the reduced 
nuclear spin degeneracy and K-level degeneracy. E u is the energy 
of the upper level and Q is the partition function of the molecule 
(with excitation temperature, T ex ). Values for the constants were 
obtained from Rohlfs & Wilson (2000) and the Jet Propulsion 
Laboratory (JPL) database (Pickett et al. 1998). To determine the 
H2 column density, we assume a 12 CO to H2 abundance ratio of 
~ 8.5 x 10~ 6 (appropriate for diffuse clouds e.g. Liszt 2007; for 
comparison, the typical values quoted for dense clouds is8xl0 -5 ). 
The peak column densities of the molecular gas iV(H2) seen to- 
wards Kepler's remnant estimated from the 12 CO emission are 3.3, 
4.4 and 9.1 xlO 19 cm -2 . 

The upper limit on the column density from the 13 CO emis- 
sion can be estimated using a similar approach to Eq. [3] with an 
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assumed linewidth and excitation temperature for the CO line set 
from the observed values for the 12 CO emission (e.g. Thompson & 
Macdonald 2003). Assuming a conversion ratio between 13 CO to 
H2 of rsj 2 x 10 -7 (e.g. Langer & Penzias 1990), we estimate the 
3 a limit on the molecular column density N(H.2) <10 20 cm -2 at 
the location where we measure both 12 CO and 13 CO (i.e. the peak 
emission in 12 CO). 

In order to confirm the validity of our LTE analysis we also 
modeled the 12 CO line emission with a non-LTE radiative transfer 
code (van der Tak et al. 2007). The code assumes an isothermal, 
homogenous medium which fills the telescope beam and solves the 
radiative transfer equation using the escape probability formula- 
tion. We modeled the integrated intensity of the 12 CO emission for 
each of the three positions given in Table [TJ assuming a kinetic 
temperature Tk=15 K and a range of different H2 densities. The re- 
sults of our modelling fully confirm our LTE analysis, confirming 
that the 12 CO emission is both optically thin and in LTE. Column 
densities derived by both methods agree within a factor of 3. The 
radiative transfer models also suggest an upper limit for the H2 col- 
umn density of ^10 20 cm -2 . 



2.2.3 Hi observations: searching for atomic structures 

To search for atomic foreground features towards Kepler, we have 
revisited the Hi observations reported in RG99, made between 
1996 and 1997, with the VLA in its CnB, C and D configurations. 
The number of spectral channels was 127, centred at Okms -1 
(LSR) with a resolution of 1 .3 km s - 1 . We studied a region of about 
25 arcmin x 25 arcmin centred on Kepler's SNR. The raw data 
were downloaded from the VLA public archive and calibrated with 
the MPS package, following standard procedures. The continuum 
was subtracted in the visibility plane (recall Fig. [2]) and the Fourier 
transform and cleaning of the images was carried out using the 
MIRIAD data processing package (Sault, Teuben & Wright 1995). 
A first data cube was constructed using natural weighting. The reso- 
lution and sensitivity achieved were slightly improved with respect 
to those reported in RG99 but the sidelobe level was unacceptable, 
so a second cube was constructed using uniform weighting which 
drastically reduced the sidelobes and improved the resolution but 
degraded the sensitivity. 

We constructed an opacity (r) cube, defined as: 

— Kt) 

where Tb (v) is the brightness temperature at velocity v and T c is 
the continuum emission. To construct the continuum image (Fig. [2]), 
the line-free channels used to determine the spectral baseline were 
combined in the u,v plane and inverted using uniform weight- 
ing. The resulting beam is 9.2arcsec x 5.5 arcsec, PA=51.6°, and 
sensitivity, 0.6 mJy beam -1 . As originally noted in RG99, the Hi 
emission is uniform in structure, this is confirmed in the opacity 
cubes and we find no evidence for a correlation of H I structures 
with the submm emission. Several trials were made to improve 
the sensitivity by Hanning smoothing the spectral profiles and con- 
volving the images to larger beams, with little success. Addition of 
single dish data from the LAB survey (Hartmann & Burton 1997; 
Kalberla et al. 2005) to the images did not add evidence of a corre- 
lation. 

Using the VLA H I and single-dish data we can estimate the 
column density of atomic gas at the location of the submm clumps: 



N = 1.82 x 10 1 * J Tdv (5) 

where v is the velocity (in kms -1 ). Integrating over the veloc- 
ity range in which large-scale Hi features are seen (—10 < 
v < 22 kms -1 ), we detect small variations in emission on the 
Hi map, with minimum and maximum column densities of 0.6- 
0.9 x 10 21 cm -2 . The atomic gas density is at least an order of mag- 
nitude higher than the molecular density estimated from the CO. 
The chopping procedure employed in the submm observations of 
Kepler's SNR can not produce the submm structures from chop- 
ping on-off these small variations, indeed chopping would have re- 
moved any emission associated with the H I as was seen in Cas A 
(Wilson & Bartla 2005). 



3 COMPARING THE GAS AND DUST TOWARDS 
KEPLER 

The high resolution 12 CO (J=2-l) map integrated over the ve- 
locity range —200 < v < 200 kms -1 , is shown in Fig. [8] with 
the outer contour of the radio image overlaid. The map has been 
smoothed with a 21 -arcsec Gaussian. Cyan crosses indicate the 
submm dust clumps, A-J; red crosses indicate CO peaks. Spectra 
at each of these locations were extracted; where CO was detected, 
the lines were fitted with Gaussian profiles with the derived prop- 
erties given in Table |2] Also listed in Tableware the fluxes in the 
submm clumps, A-J, determined by placing an aperture of radius 
23 arcsec on the synchrotron- subtracted 850-/xm image. The sum 
of the fluxes in these apertures is 0.38 Jy (cf. 0.58 Jy in M03). The 
submm clumps outside the remnant are located on average, in re- 
gions of higher CO emission than those inside the remnant, how- 
ever a significant fraction of the submm flux comes from clumps 
where there is no, or little CO emission. Comparing the distribution 
of the CO and submm emission spatially, we can see that molecu- 
lar clouds are observed at the location of submm clumps G and H 
(combined flux of 94 mJy) and fainter emission at clumps A, C, I 
and J (combined submm flux 135 mJy). We do not detect CO emis- 
sion at the locations of the submm clumps B, D, E and F (combined 
submm flux 153 mJy). 

The CO column densities towards these clumps is a better in- 
dicator of the quantity of line-of-sight emission. At apertures A- 
J, these range from <10 14 -10 15 cm -2 with molecular gas densi- 
ties 10 19 -10 20 cm -2 . The gas column densities estimated from the 
submm fluxes range from 10 21 to 10 22 cm (with gas-to-dust ra- 
tio of 160); as confirmed by the column densities estimated from 
the (lower-resolution) 100/xm IRIS maps (Miville-Deschenes & La- 
gache 2005). This a factor of 10-100 times higher than that esti- 
mated from the CO data. If the submm clumps were related to the 
same molecular gas as traced by the foreground CO emission, the 
conversion ratio (H2)/ CO in these clumps would have to be less 
than 10 - 7 which is far lower than the values (>10 -5 ) expected for 
clouds with densities of 10 21 -10 22 cm -2 (e.g. Listz 2007). Con- 
versely, using the appropriate conversion factor for the molecular 
gas densities estimated from the submm clumps, we would expect 
to observe CO column densities (CO) >10 16 -10 17 cm -2 at the 
location of the submm emission. This suggests that the molecular 
gas contribution to the submm emission towards Kepler is negli- 
gible, even where CO is co- spatial with the submm clumps on the 
sky. 

Are there any physical effects that could result in our underes- 
timating the column density of the gas? In dense molecular cores, 
CO and other molecules may become depleted from the gas phase 
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Figure 8. Integrated 12 CO(J=2-l) emission towards Kepler's SNR over the velocity range —200-200 km s _1 smoothed by a 3kms _1 Gaussian. The crosses 
indicate regions where spectra were taken across the CO clouds and the remnant. Cyan crosses show the locations of the submm clumps A-F, red crosses are 
positions chosen with strong CO emission for comparison. The integrated spectra from apertures at the locations of the submm clumps and at the CO clouds 
are shown alongside over the velocity range —50 < v < 50kms _1 (unsmoothed). The greyscale ranges from —1.4 to 1.9 Kkms -1 . 



by freezing out onto the surfaces of dust grains (e.g. Redman et al. 
2002). The timescale for molecular freezeout is tightly correlated to 
the gas density and only becomes significant for H 2 densities >10 5 
cm -3 (Caselli et al. 1999). If this were the case then we would ex- 
pect to see evidence for higher gas densities in the form of higher 
optical depths and the presence of rarer isotopologues of CO. The 



properties of the CO emission that we have detected are strongly 
suggestive of diffuse clouds in which molecular depletion does not 
play a significant role in controlling the abundance of CO. 

In diffuse molecular clouds, the CO abundance depends on the 
overall H 2 column density (Liszt et al 2007), likely due to the in- 
creased UV photodissociation rates for CO in low density gas (van 
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Name Coords T peak v c fwhm r.m.s. /(CO) 7V( 12 CO) iV(H 2 ) S 8 50 

AR.A. A Dec. K (kms" 1 ) (kms- 1 ) K (Kkms" 1 ) (10 14 cm- 2 ) (10 19 cm" 2 ) (mJy) 



COS 


+221.9 


+95.9 


0.63 


4.06 


1.14 


0.05 


1.12 


6.0 


7.1 










0.51 


11.17 


0.46 


0.05 


0.25 


1.1 


1.3 




COS 


+146.3 


+82.7 


0.67 


4.05 


1.10 


0.05 


0.78 


4.3 


4.9 




COS 


+51.1 


+192.5 


0.46 


4.30 


1.00 


0.06 


0.49 


2.7 


3.1 




CO N 


+96.9 


-190.9 


0.34 


11.48 


1.17 


0.05 


0.42 


2.2 


2.7 




CON 


+47.6 


-153.5 


1.14 


11.12 


0.18 


0.06 


0.81 


4.4 


5.1 




A 


-83.3 


+23.3 


0.14 


11.4 


0.96 


0.05 


0.16 


8.6 


10.2 


28.8 ± 3.5 


B 


-17.5 


+71.8 








0.05 


<0.15 


<0.8 


<1.0 


23.0 ± 2.9 


C 


+69.3 


+110.3 


0.13 


11.21 


1.59 


0.05 


0.22 


1.2 


1.5 


41.5 ±4.6 


D 


+101.5 


+93.3 








0.05 


<0.15 


<0.8 


<1.0 


45.8 ±4.9 


E 


+109.9 


+47.3 








0.05 


<0.15 


<0.8 


<1.0 


62.1 ±6.6 


F 


+94.5 


+58.3 








0.06 


<0.18 


<1.0 


<1.1 


21.8 ±2.7 


G 


+65.1 


-57.0 


0.20 


11.27 


0.87 


0.05 


0.19 


10.2 


1.2 


30.6 ± 3.6 


H 


+70.7 


-104.7 


0.55 


-1.59 


0.48 


0.05 


0.28 


1.6 


1.8 


63.8 ± 6.7 








0.55 


11.36 


0.76 


0.05 


0.42 


2.2 


2.7 




I 


+21.7 


-71.7 


0.23 


11.71 


0.68 


0.05 


0.20 


1.1 


1.2 


31.4 ±3.3 


J 


-7.7 


-118.7 


0.27 


11.70 


0.72 


0.05 


0.21 


1.1 


1.3 


33.5 ±3.8 



Table 2. Properties of the detected 12 CO ( J=2-l) line emission towards Kepler's SNR. Columns show (1) name (2) offset coordinates from centre in arcsec 
(RA = 17 h 30 m 41 s .3 and Dec. = -21°29 / 29 // where postive R.A. and Dec. indicates a move to the East and South); (3) peak T MB (4) central velocity 
(5) velocity width (fwhm) (6) integrated intensity (7) rms (8) and (9) column densities of CO and H2 gas estimated using Eq.[3]and (10) the submm flux 
measured in the apertures A-J after synchrotron subtraction. Errors are estimated using the procedure outlined in Dunne & Eales (2001). 



Dishoeck & Black 1988). However, in order to match the column 
densities indicated by the submm emission, the density of the gas 
would increase to a point where photodissociation is largely unim- 
portant except at the exterior of the cloud. The CO to H 2 abundance 
ratio for gas with an H2 column density of 10 21 -10 22 is close to that 
expected for a dark molecular cloud and not consistent with the col- 
umn densities of CO that we observe. We thus conclude that neither 
depletion or photodissociation are likely to have a significant effect 
upon our derived column densities. 



3.1 The revised dust mass in Kepler's remnant 

The IR-submm spectral energy distribution (SED) of Kepler's 
SNR, using the synchrotron- subtracted submm fluxes is shown in 
Fig. [9] We fit a two-component modified blackbody to the SED 
which is the sum of two modified Planck functions, each with a 
characteristic temperature, T w and T c (Eq.[6]): 

S„ = N w x v^B(v, T w ) + N c x v*B{y, T c ) (6) 

where N w and N c represent the relative masses in the warm 
and cold component, B(v,T) is the Planck function and /3 is the 
dust emissivity index. The model was fitted to the SED (constrained 
by the 12-/xm flux) and the resulting parameters (N w , N C ,T W ,T C , 
/3) which gave the minimal x 2 were found. 

The solid curve shows the fit to the SED with 5ioo = 5.6 ± 
1.4 Jy (the average from Arendt 1989 and Saken, Fesen & Shull 
1992) whilst the dot-dashed curve fits the SED with the lower 
100-/im flux of 2.9 Jy (estimated from pointed observations, Braun 
1987). The higher 100-/xm flux predicts the lower dust mass for 
the SED model. Both fits include the recent Spitzer fluxes at 24 
and 70 /xm (Blair et al. 2007) but the 160/im upper limit is not 
included in the fit due to uncertainties with the background at this 
wavelength. Neither fit to the IR data alone would provide evidence 
for a cold dust component, highlighting the importance of accurate 
submm fluxes when estimating the dust mass. Note that the SED is 
equally well fit in both cases and to distinguish between the two, 




Wavelength (//.m) 

Figure 9. IR-submm SED for Kepler's SNR, comparing the best fits with 
two different 100-/mi fluxes, including Spitzer data and submm points (this 
work). The solid line represents the best fit with S100 ~ 5.9 ± 1.4 Jy. The 
dot-dashed curve represents the fit to the lower 100-/im flux, 2.9 ± 0.9 Jy. 
The upper limit shown by the arrow is the total integrated flux obtained 
from Spitzer at 160 pm (Blair et al. 2007). 



more accurate fluxes are needed around the peak of the cold emis- 
sion at 200-400 pm. 

In addition to using the range of IRAS 100-/im fluxes in the 
literature, we also applied a bootstrap analysis to our SED fitting to 
determine the errors on the model. The photometry measurements 
were perturbed according to the errors at each wavelength (using 
the higher 100-/im flux which produces the lowest mass). The new 
fluxes were fitted in the same way, the SED parameters and dust 
masses recorded, with the procedure repeated 3,000 times. The me- 
dian results along with errors (estimated from the 68 per cent con- 
fidence intervals) are listed in Table [3] 

The dust mass is calculated from the SED model using Eq. [7] 
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Bootstrap parameters Dust mass (Mq) 

T W (K) T C (K) Aci^kg" 1 ) ACaCm^g" 1 ) 

2.0 ±0.2 88^ 16 ±1 0.12+°;° 2 ^toio 



with the remnant ranges from 0. 1-1.2 Mq, depending on the ab- 
sorption coefficient. This value is 100 times larger than seen by 
Spitzer and concurs with the results for Cas A (Dunne et al. 2009) 
that supernovae are significant sources of dust. 



Table 3. The best-fit parameters for Kepler's SNR estimated from the boot- 
strap technique. The dust mass is the median of the distribution from the 
bootstrap technique with the error quoted from the 68 per cent confidence 
interval, ki = 0.7 m 2 kg -1 (appropriate for SNR dust, Dunne et al. 2009) 
and K2 = 0.07 m 2 kg -1 (typical ISM values) which gives the maximum 
mass. 



APPENDIX A: 
CHOPPING? 



SPURIOUS STRUCTURE DUE TO 



M d 



^850 



N w 



B(S50 fim,T w ) B(850/im,T c ) 



(7) 



where k,85o is the dust absorption coefficient and the distance, D 
is taken to be the lower published limit of 3.9 kpc (Sankrit et al. 
2008). To determine the dust mass, we take two extreme values for 
k\ (i) ^sso = 0.07 m 2 kg -1 , typical quoted value for grains in the 
diffuse interstellar medium (D03 and references therein) and (ii) 
^850 = 0.7 m 2 kg -1 , as required for the submm polarimetry obser- 
vations of Cas A (Dunne et al. 2009). This is similar to the values 
(~ 0.4 m 2 kg -1 ) predicted by the supernova dust model in Bianchi 
& Schneider (2007). The range of dust masses obtained from the 
SED model (Eqs.[6l&[7]) is therefore 0.1-1.2 M© depending on k. 
Note that k is the largest source of error when estimating the mass 
of dust from IR/submm emission. 



4 CONCLUSIONS 

We have presented information on the reduction and analysis of 
SCUBA submm data of Kepler's SNR including a more accurate 
subtraction of the synchrotron component. The residual fluxes are 
slightly lower than those previously reported. We find that the ring- 
like structure seen in the submm cannot be reproduced by chop- 
ping on-off large scale structures. Large-scale FIR/submm imaging 
is crucial to determine whether the observed submm structures are 
unique to the SNR, or are common in this region. Such observa- 
tions will soon be possible with SCUBA-2 and the Herschel Space 
Observatory. 

We investigated whether foreground molecular or atomic 
structures could be responsible for the submm emission conclud- 
ing that: 

• There are three molecular clouds in the vicinity of Kepler's 
SNR, located on the periphery of the remnant in the north and 
south-east and extending further out in the south. The clouds are 
faint and cold (TJ <1 K, Av < 1 kms -1 ). Molecular gas column 
densities towards the remnant are estimated to be <10 20 cm -2 . 

• 13 CO line emission is not detected at the locations of the peak 
12 CO emission. The 3-cr upper limit on the molecular column den- 
sity is <10 20 cm -2 indicating that the clouds are diffuse and opti- 
cally thin. This is confirmed with an independent analysis. 

• The column densities estimated from the submm are 10-100 
times higher than the column densities estimated from the molec- 
ular gas. This difference cannot be explained by reasonable varia- 
tions in the conversion factor between CO and H2 emission nor can 
it be explained by depletion or photodissociation effects. 

The molecular contribution to the submm emission towards Ke- 
pler's remnant is therefore negligible. The dust mass associated 



The observing technique used when obtaining the original SCUBA 
data ( §2.11 could affect the apparent distribution of submm emis- 
sion around Kepler's SNR if we had chopped onto structure un- 
related to the remnant. We modeled the effect of chopping onto 
interstellar structures not associated with the remnant using the 
chop positions and throws chosen for Kepler (Fig. lAll) . The sim- 
ulations were performed by sampling the sky model at the appro- 
priate on-off source chop locations, white noise was added before 
reconstructing the on- source flux measurements from the three dif- 
ferent chop/nod positions. We included different input models of 
the sky, here we discuss the results using (i) large scale interstellar 
submm structures (with the CO emission as a guide, Fig. lA2a ) and 
(ii) large scale interstellar submm structures including a SNR ring- 
like structure at the location of Kepler (based on the radio image, 
Fig.|A2H). 

The results are shown in Figs. IA2b and \A2b along with the 
original SCUBA data in the same colour scale (c and f). We see 
some spurious structure in the simulated output for case (i) as a 
result of chopping on to the bright extended structure in the south 
seen in the input sky. The effect of this chop is diminished when 
averaged over the two off-source chop positions (jiggles three and 
five). The strong ring-like emission structure seen in the SCUBA 
data (c) is not reproduced in the simulation, suggesting that the 
ring cannot be a result of chopping on or off large scale background 
structures. This is supported by the output from case (ii) which in- 
cludes a SNR model in the input (e). The ring-like structure sur- 
vives after chopping on-off the larger surrounding structures. We 
conclude that the emission in the SCUBA map is due to a ring-like 
source of submm radiation and is not an artifact of the observing 
technique. 
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